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Abstract 

Some implications of the proposal that flavor nondiagonal couplings of neu- 
trinos to gravity might resolve the solar neutrino problem are considered in 
the context of three neutrino flavours. Using an adiabatic approximation 
iso-SNU curves are calculated for the neutrinos which most likely contribute 
to the observed solar neutrino deflciency. We show that results obtained with 
two-flavor models can be recovered, and discuss the effects of the addition of 
a third flavor. Some results are obtained for this case, and these values are 
compared with recent experimental data. 
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1 Introduction 



The Solar Neutrino Problem (SNP) is a mystery which mars our current 
understanding of stellar astrophysics. There are now four experiments [|l], ^ 
1^, ^ employing different detection techniques that consistently find a deficit 
in the observed z/g flux relative to that predicted by various solar models 
(see e.g. those listed in 0). The origin of this discrepancy is unknown, 
but various possible solutions based on particle physics techniques have been 
proposed. 

One possible class of solutions to this dilemma invokes oscillations be- 
tween the various flavors or generations of neutrinos. Perhaps the most 
elegant solution of this type proposed to date is the Mikheyev-Smirnov- 
Wolfenstein (MSW) Mechanism ^; here, the Hamiltonian becomes off- 
diagonal under a change of eigenbasis for the neutrinos, resulting in flavor- 
oscillations, much akin to those found in the quark sector. However, the 
MSW effect requires the neutrinos to possess a mass eigenbasis, and hence 
implies massive neutrinos. 

In 1988, Gasperini conjectured that a similar oscillation phenomenon 
could take place, while still ensuring massless neutrinos. In this mechanism 
each neutrino flavor is assumed to couple differently to gravity, permitting 
flavor oscillations at the expense of violating the equivalence principle. This 
mechanism has been investigated in a number of papers (see e.g. |P, ^ p!T|), 
and has recently been referred to as the VEP mechanism [|12] . The results of 
these studies indicate that, although strongly constrained by existing data, 
the VEP mechanism is at present a phenomenologically viable means of gen- 
erating neutrino oscillations. 

All investigations of the VEP mechanism to date have restricted consider- 
ation to 2-flavor oscillations. However, at our current level of understanding, 
we know there to be at least three generations of neutrinos. While two-flavor 
oscillation models are elegant and mathematically simple, it remains a mys- 
tery whether or not the r-neutrino {vr) plays a crucial role in determining the 
observed v counting rates and fluxes. Several studies 0, |15|, [16|, have 
shown that, for certain three-flavor analysis can and does contribute 

interesting physics to the SNP. 

In this paper we consider an adiabatic extension of the VEP mechanism 
to the realistic three-flavor case. The width of the parameter space region 
found in [jl2| is reproduced for a two- flavor limit, and the size of the region is 
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found to shift and/or widen for certain values of the Ve mixing angle. 

Spectra for the solar neutrinos are produced for differing values of the 
mixing angles (for fixed violation parameters A/21 and A/31), both for a fixed 
neutrino event rate, as well as for a fixed event rate as would be detected 
by the ^''C/ experiment [i.e. combined counting rate for ^B + '^Be neutrinos). 
These are compared with the predicted flux in the two-flavor limit. 



2 Review of Three-Flavor MSW Neutrino 
Oscillations 

Before proceeding to the gravitational case, we review the MSW mechanism 
for three flavors. Consider the neutrino eigenbases 



Ww) 



I > V'w) = R3 Wm) , (1) 



where {i^w) is the weak eigenbasis and |z/m) the mass eigenbasis. The matrix 
-R3 is the 3x3 leptonic analogue of the CKM mixing matrix, which we 
parametrize as 

and is proportional the the identity matrix in the case of massless neutrinos. 
Here ip, 0, cu are three mixing angles, and the SU{3) generators. For 
simplicity we have taken R-^ to be real, which implies that there are no CP 
violations in this sector (i.e. the CP- violating phase 6 = 0). Including 
matter effects, and (without loss of generality) assuming ma > m2 > mi, the 
Hamiltonian in the mass eigenbasis becomes 

H' = 

^ ml + ml + Acl-AN As2^ ^{l-P)j2 \ 

\ m^ + m^ + A4 + A jv A 52^^(1 + /9)/2 

^ As2^^{\-^)I2 As^^{l + l3)/2 2ml + 2Asl 

■ (3) 
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where = cos0, etc. We follow the definitions of [|T^] with 



= Ucl-A,C2.f + Alsl]'/' (4) 
/5 = -^{Acl-A,c,^r (5) 

where = ml — ml, and A = A{r) is the matter interaction term, 

A{r) = 2V2GFNe{r)E . (6) 

which describes the interaction of i/e-neutrinos with solar matter. 

From this point on, we make the following substitutions for the mixing 
angles of (0): u = 9i2, (p = O13, ip = 6*23. Assuming that the the state 
transitions are adiabatic, we can approximate the z/e survival probability 

P,^^,^(r,E) to be 

Pu^^uAr,E) = {1-Plz2)0{A~A2i) + Plz2{1-Plz3)0{A2i-A) 

xQ{A- A31) + PlziPlz3 e(A3i - A) . (7) 

Here Q{x) is a step function I, with 

@C^^ ^ 1 1 Vx < 

1 otherwise ' 

and PLZi is the Laudau-Zeener "jump" probability for the z/g z/j transition, 
with 

/ cos 29ii \ ^ TT 2 n \ / \ 

tti = 2TXKi . ^ , A = 7T'^i(l - tan On) . (10) 

y sm Zvii J / 

For z/g transitions, i = 2, and i = 3 for z/g ^ z^r- 



^The step function is a good approximation for the sines and cosines of the matter- 
enhanced mixing angles, in the case of smaU vacuum angles The results obtained 
using this form of the survival probability are in excellent agreement with the full form of 
the (adiabatic) probability. 



3 



3 Gravitationally-Induced Neutrino Oscilla- 
tions 



We now modify the equations of the previous section to model neutrino 
oscillations induced by the non-zero difference Afij between the gravitational 
coupling of the different flavors. Consider first the 2- flavor scenario. Let {i^w), 
|z/g) be the respective weak and gravitational eigenbases 

where the two are related by an SO (2) rotation R = R{9g) (the "mixing" 
matrix) , 

\vw) = R\i^g). (12) 
The gravitational eigenstates obey a variation of the Dirac Equation 

[i{e/i^d,)-Cm\^G) = ^. (13) 

where C{'f) is some function of the gravitational potential, which arises due 
to the non-zero violation of the equivalence principle. For a spherically sym- 
metric metric (choosing angular coordinates 6' = = as the trajectory of 
the neutrino), the vierbeins in (p!3D reduce the evolution equation to fl^ 



= H \vg) . (14) 

where 

H, = -2\<P{r)\E,{l + f\) . (15) 
are the diagonal components of the Hamiltonian in the gravitational eigen- 



basis [0. Here, <C 1 is a (flavor-dependent) parameter which determines 
the magnitude of the violation of the Einstein Equivalence Principle (EEP). 
Note that for (i.e. first generation matter), / = 0, which differs from the 
notation of [jl2| (in which / = 1 represents first-generation neutrinos, which 
is equivalent to saying that the term (1 -|- /) = 1 in (0)). Under the change 
of basis \vg) \i'w)i the Hamiltonian becomes off-diagonal, 

= R-^HR\vw) , (16) 

dr 
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and with the inclusion of weak interactions in the v^, sector, we obtain the 
effective Hamiltonian 



E' 



R-^HR 



We 



(17) 



^GFNJr)-E 



E\ 



r)\Afcos2eG E\ 
r)|A/sin2^G 



(r)|A/sin2^G 




These matter-enhanced oscillations are resonant for complete flavor con- 
version, which implies the energy condition 



E 



V2GFNe{r) 
2|0(r)|A/cos2^G 



(18) 



where Gp is the Fermi Coupling Constant, Ne{r) the radial electron number- 



density, and 0(r) the dimensionless (solar) gravitational potential. 
A/ = /2 — /i. In addition, we require the adiabaticity condition 



V2Gf {Ne)res tan22^G 





1 dNe _ ld4 






Ne dr (p dr 


res 



> 1 . 



Also, 



(19) 



Iso-SNU curves [i.e. curves of constant neutrino-counting rates) for 2- 
fiavor gravitationally-induced neutrino oscillations have been carried out in 
|JTl|, and later repeated in [1^. Two distinct regions of overlap in pa- 
rameter space were discovered for the given experimental data [0], namely 
a highly non-adiabatic region (small mixing region), and an adiabatic "large 
mixing" region. However the authors in ^ have argued that these regions 
are too small, statistically speaking, to be kept as viable solutions to the 
SNP for equivalence principle violating neutrino flavors. 



Extension of VEP Mechanism to Three 
Flavors 



Turning now to the three flavor case, we substitute 



1 

2E 



( 















mq 



2E|0(r)| 




(20) 
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which amounts to replacing the mass-squared difference in the probabihties 
(1^ with an EEP violation term, 



where i,j G [1,3], i ^ j. For the survival probability, just the 

2 = 2, 3, j = 1 terms come into play, and we need only concern ourselves 
with the mixing angles Qyi-, On- 

Of course the entire parameter space (A/21, A/31, dyi, Qn) is not fair game 
for this analysis. In particular, we must limit the space to only the areas in- 
volving resonances, and which are allowed by adiabatic processes. The phys- 
ical boundaries on the A/ parameter are easily calculable, and have in the 
two-flavour case been found to be |]12[ A/ G (10~^^, 10~^^). The same range 
holds for the parameters A/jj. This implies that a v^. with |-EA/| > 10~^^ 
MeV will not undergo resonance, and will emerge from the sun unscathed, 
retaining its identity The other A/ boundary is similar, but for the 

lower resonance. 

This determination of resonant behavior provides us with valuable infor- 
mation. There are now three distinct regions of parameter space which are 
of physical interest: 




Am?- 



(21) 



1. A/31 > 10-12, A/21 G (10-18,10-12): 



— Here, we should expect large dependence on — »• z/g transi- 
tions [i.e. A/21, ^12), as the A/31 parameter puts all transi- 
tions to i^T- above the upper resonance. In the limit ^13 — *• 0, we 
should expect to recover the two-flavor EEP violation model. 



2. A/3le((10-l^lO-l2),A/,, <io- 



— This range represents the opposite of the previous one, where 
the Vf. z/^ transitions are not resonant. Hence, we expect to 
see a two-flavor limit with strong dependence on the v^. Vt 
transition parameters (^13, A/31). 



3. 10-18 < A/21 < A/31 < 10- 



— The most interesting behavior is expected to occur here, where 
both parameters are sandwiched into the allowed A/ param- 



eter range. Similar to the mass MSW mechanism we 
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should expect to see dependence on all four parameters in 
question. 

The following section summarizes our analysis of the physically relevant 
cases #1 and #3. 



5 Results 

Using the neutrino spectra and solar mass data from |p, as well as the ^^Cl 
and ''^a cross section data from [1^, we have computed iso-SNU rates for 



case #1 as listed above. The counting rate for a given detector material 
X with energy threshold Emm and absorption cross-section a^i^E) is 

a=/ drr\{r) dE <P{E)a^{E)P,^^,Xr.E) , (22) 

for neutrinos with energy spectrum and maximum energy E^ax- The 

function ^(r) represents the fraction of neutrinos of a given type produced at 
radius r. The rates obtained by numerical intergration of (^21) are essentially 
commensurate with those obtained in ref.|p and others {e.g. see |T2|) for the 
"large mixing region" (to within a few percent). The neutrinos which were 
assumed to contribute the largest counting rates for each type of detector 
include, in order of abundance: 



Note a slight discrepancy in the upper limit of the region as given in [|T2 
It is believed that any discrepancies could well be the result of numerical 
errors during calculation, or from minute variations in data for the specific 
neutrinos. We show in Table 1 the (Sex) adiabatic edges of the parameter 
space in question, according to the rates as given in [O: 



Cnome. = 2.55 ± 0.25 SNU 
Cgallex = 79 ±11.7 SNU 
CsAGE = 73 ± 19.3 SNU . 
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The edges of the parameter space region for the Homestake ^^C/ detector 
are roughly equal to the bounds in question, while the 3a range for the 
''^a experiments is much larger. Hence, the overlap range is determined 
predominantly by the ^^C/ counting rates, which have been calculated to be 
sin^26'i2 G (0.60, 0.89). This is found to be in close agreement to the range 
width obtained by p[ of sin2 20i2 G (0.60, 0.90). 

As previously mentioned, the most interesting physics is expected to 
emerge from case 7^3, where both resonances become crucial in determi- 
nation of the neutrino behavior. We have obtained some preliminary results 
for various values of the parameters, using the data for the ^^C/ detector. 
We find that upon inclusion of the third flavor, certain parameters may be 
chosen such that the "large mixing region" can be widened. To see the mag- 
nitude of the widening, we compare the width of the region for 2 flavors to 
the width using 3 flavors (with two choices of ^13). For the three-flavor model 
{i.e. Ovi ^ 0), we choose A/31 = 10-^^ 

Table 1: Width of Allowed Parameter Space for Large Mixing Region as a 
function of sin^ 6*13 



sin^ 6^13 


Width (sin^ 2^12) 





0.60 - 0.89 


0.01 


0.65 - 0.94 


0.05 


0.64 - 0.93 


0.10 


0.63 - 0.93 


0.40 


0.57- 0.91 



From table 1 we see that the known parameter space overlap for the Ve 
i>fj_ transitions is shifted and expanded somewhat with the addition of the Ur. 
Further study of the adiabatic behavior of the neutrinos with respect to the 
SNP may yield a nice set of parameters which maximize the allowed region. 
Also, we may narrow down our possible choices of violation parameters and 
mixing angles through a comparison of detected to theoretically calculated 
neutrino spectra. Fig. 1 shows the calculated flux for two of the parameter 
sets listed in Table 1, for a fixed SNU rate of neutrinos (i.e. C = 1.38 
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SNU). Note the subtle shifts in shape of each reduced spectrum for the lower 
half of the energy range, i.e. below approximately 6 MeV (Fig.l), as well as 
the flattening of the flux curve for increasing ^13 in Fig. 2. These differences 
could possibly be differentiated by present-day or near-future detectors (see 



e.g. some discussions in |]T3|, |T6|, p!8|). 

A VEP analysis incorporating all three generations of neutrinos seems to 
yield a greater mobility in parameter space. Such an analysis is quite timely, 
due to the advent of such new ventures as the Sudbury Neutrino Observatory, 
for more precise measurements of neutrino fluxes can help determine the right 
set of parameters. As we can see from Figs. 1 and 2, there is greater range 
of parameters from which to choose such that we can exactly match a flux 
curve. Since SNO will be able to detect the full solar neutrino flux (flavor- 
independent via neutral current interactions), studies of such spectra could 
also help determine the reduction mechanism at work, and hence possibly 
distinguish between a VEP- or MSW-based phenomenon. Also, an analysis 
of the non-adiabatic transitions is currently underway, to gauge what effect 
this will have on the full parameter space range. Since the small mixing 
region is dominated by highly non-adiabatic transitions, it is believed that 
this could represent a larger portion of the expanded three-flavor space. 
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Figure Captions 



Fig. 1: Plot of reduced neutrino fluxes {cm~'^s~^MeV~^) v.s. neu- 
trino energy (MeV) for flxed counting rate Csb = 1-38 SNU showing vari- 
ation of spectrum shape with choice of parameters. Upper (large) flux is 



unreduced SSM B flux from |jT9[, while reduced flux curves for 3-flavor 
model are calculated with A/21 = 10~^^, A/31 = 10~^^, and 

• sin'^2eu = 0.685 ; sm^^ig = 0.40 (dotted) 

• sm22^i2 = 0.70 ; sm^^is = 0.05 (sohd) 

These are compared with predicted 2-flavor flux (A/31 > 10~^^), with 

• sin'^2eu = 0.697 ; sm^^ig = (dashed) 



Fig. 2: Plot of reduced ^B neutrino fluxes v.s. neutrino energy for flxed 
counting rate as would be detected by ^'^Cl experiments (i.e. ^B + "^Be), 
Cci = 2.00 SNU. Reduced flux curves for 3-flavor model are calculated with 
A/21 = 10-^^ A/31 = 10-l^ and 

• 3171^2012 = 0.63 ; sm2^i3 = 0.40 (dotted) 

• sin^29i2 = 0.693 ; sm^^ig = 0.05 (solid) 

These are compared with predicted 2-flavor flux (A/31 > 10~^^), with 

• sm22^i2 = 0.701 ; sin^^is = (dashed) 
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Fig.l: B8 flux; Fixed Counting Rate 




This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9501340vl 



Fig.2: B8 flux; Fixed CI Counting Rate 




This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9501340vl 



